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Recent successes in manufacturing of atomi- 
cally thin graphite samples [1] (graphene) have 
stimulated intense experimental and theoretical 
activity [2, 3]. The key feature of graphene is the 
massless Dirac type of low-energy electron exci- 
tations. This gives rise to a number of unusual 
physical properties of this system distinguish- 
ing it from conventional two-dimensional met- 
als. One of the most remarkable properties of 
graphene is the anomalous quantum Hall effect 
0) la 0> B @|« It is extremely sensitive to the 
structure of the system; in particular, it clearly 
distinguishes single- and double-layer samples. In 
spite of the impressive experimental progress, the 
theory of quantum Hall effect in graphene has not 
been established. This theory is a subject of the 
present paper. We demonstrate that the Landau 
level structure by itself is not sufficient to deter- 
mine the form of the quantum Hall effect. The 
Hall quantization is due to Anderson localization 
which, in graphene, is very peculiar and depends 
strongly on the character of disorder [9] . It is only 
a special symmetry of disorder that may give rise 
to anomalous quantum Hall effects in graphene. 
We analyze the symmetries of disordered single- 
and double-layer graphene in magnetic field and 
identify the conditions for anomalous Hall quan- 
tization. 

As was discovered 10] in 1980, the Hall conductivity 
a xy of a 2D electron gas in a strong transverse magnetic 
field develops plateaus at values quantized in units of 
e 2 /h. This phenomenon is the famous integer quantum 



Hall effect (QHE) [ll| — one of the most fascinating 
quantum effects in the condensed matter physics. 

The experimentally measured Hall conductivity of 
single-layer graphene [1, IH, 0] is quantized taking the odd 
multiples of the quantum 2e 2 /h (here the factor of two 
is due to the spin degeneracy) 

a xy = (2k + l)2e 2 /h, 



k € 



(1) 



In double-layer samples, the quantum Hall plateaus oc- 
cur at even multiples of 2e 2 /h excluding k = 0. Ow- 
ing to this unusual quantization, the Hall measurements 
are widely used in modern experiments for characteriz- 
ing the graphene samples. Remarkably, the signatures of 
Hall conductivity quantization in graphene were recently 
observed even at room temperature [7( . 

A simple argument in favour of the odd QHE ([I]) in 



a single graphene layer 12j is based on the structure 
of Landau levels for two-dimensional massless electrons 
[ijj]. In clean graphene, the energies of Landau levels 
are ejv = huj c sgnN y/\N\ with lu c = vo\/2eB/hc and 
N € Z. The plateaus of the Hall conductivity are then 
identified with its classical values a xy — n e ec/ B at con- 
centrations n e corresponding to an integer filling factor 
ri = n e hc/eB, that is, to an integer number of filled Lan- 
dau levels. This consideration is further extended by the 
calculation of the Hall conductivity in the presence of 
disorder within Boltzmann [l2j] or self-consistent Born 
approximation (l4| . i.e., for disorder-broadened Landau 
levels. 

However, the spectral gaps in the density of states be- 
tween separated Landau levels do not lead to the QHE. 
Indeed, while the dependence of the Hall conductivity on 
the Fermi energy is quantized in the clean system, this 
is not the true QHE. The point is that the Fermi level 
itself is not a smooth function of the density: it jumps 
between the fully occupied and empty Landau levels with 
increasing density. As a result, the density dependence 
of <j X y (which is measured in experiments) shows up no 
steps and no plateaus, i.e. no QHE. Including the Landau 
level broadening by disorder leads only to the magneto- 
oscillations of cr X y(n e ) but not to its quantization [l~5| . 
In fact, (i) the QHE does not require quantization of 
the density of states at all; (ii) the position of the QHE 
plateaus transition does not necessarily correspond to the 
center of Landau level; (iii) the crucial ingredient respon- 
sible for the Hall quantization is the disorder-induced An- 
derson localization (for review, see [HI, [H, 17|). 

The existence of the odd-integer ((T|) QHE in graphene 
also requires a more rigorous justification in view of quan- 
tum interference effects that are essential in any two di- 
mensional system, including graphene. Once disorder is 
fully taken into account, the quantization of Hall con- 
ductivity is exact and the transition between quantum 
Hall plateaus becomes a quantum phase transition with 
universal critical properties. This immediately shows the 
non-universality of the result |T]). Indeed, if the disorder 
in graphene is of a generic form and does not possess 
any special symmetry, then the Dirac nature of excita- 
tions will be completely lost at large length scales. This 
is exactly what occurs in graphene with a generic (pre- 
serving only the global time-reversal symmetry) disorder 
at zero magnetic field B. In such a system, localization 
yields vanishing conductivity with lowering temperature 
The critical properties of the generically disor- 
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dered graphene will not differ from those for any other 
two-dimensional system. The quantized Hall conductiv- 
ity will then take all integer multiples of 2e 2 /h rather 
than the odd series Eq. (TT]). Furthermore, a xy = at 
the Dirac point due to the particle-hole symmetry. The 
conventional theory of the QHE [2(| predicts complete 
localization [a xx = 0) and the plateau in the Hall con- 
ductivity instead of the quantum Hall transition in this 
situation. Therefore, the observation of the odd quan- 
tization ([1]) is a striking experimental result calling for 
theoretical explanation. 

The only reason for a non-standard quantization of the 
Hall conductivity is the presence of some special symme- 
try that is preserved by disorder and thus changes the 
critical behavior of the system. Unconventional trans- 
port and localization properties of graphene with spe- 
cial symmetries of disorder at B = were studied in 
Refs. 21, 2 p| fs ee also earlier works on disordered Dirac 



fermions ). However, the quantum localization ef- 

fects (and hence most of the peculiarities arising from the 
symmetry of disorder) were discarded in most analytical 
studies devoted to the QHE in graphene 13, 3, 25, 2t| . 
Recent numerical simulations of disordered graphene in 
magnetic field 27, 28, [29| have indeed shown that the 
result is sensitive to symmetry properties of disorder. 

In this paper we develop the theory of the integer QHE 
in graphene. We carry out the symmetry analysis and 
identify the situations when the QHE is anomalous. 

Graphene: model and symmetries. — We start 
with the effective Hamiltonian for the clean single-layer 
graphene in external magnetic field 



H 



V T 3 (T 



(2) 



Here the Pauli matrices a% and Tj operate in the space 
of two sublattices A and B and two valleys, K and K', 
of the graphene spectrum, respectively. The full symme- 
try classification for this Hamiltonian in the absence of 
magnetic field was developed in Ref . 2l| . When the mag- 
netic field is applied, time-inversion symmetry is broken 
and we are left with (i) an SU(2) isospin symmetry in 
the space of valleys [3fJ, generated by A XiV = 0-3x1,2 and 
A z = CT0T3 and (ii) an additional discrete chiral symme- 
try Co that arises exactly at zero energy: H = —(T3H0-3. 
Further, we denote C X)3/jZ the combinations of Co trans- 
formation with the isospin rotations. 

We first consider the situation when all chiral symme- 
tries are broken. This always happens when the Fermi 
energy is shifted away from the Dirac point by the gate 
voltage. At zero Fermi energy, the chiral symmetry can 
be violated by e.g. any potential disorder. In this case we 
have only two possibilities with respect to the symmetry: 
decoupled (A z preserved) or mixed (A z violated) valleys. 

Decoupled valleys: odd quantum Hall effect. — 
We start with considering the case of decoupled valleys. 
A physical realization is any disorder smooth on the scale 



of lattice spacing. The isospin of electrons (valley index) 
is hence preserved. We will show that it is the isospin 
symmetry that is responsible for the odd quantization 
Eq. ([T]) . Indeed the isospin degeneracy implies the quan- 
tization of Hall conductivity with the step 4e 2 /h (the 
factor 4 accounts for 2 degenerate spin states and 2 in- 
dependent valleys). Then to prove the validity of Eq. 
([1]), it suffices to establish the quantum Hall transition 
at zero filling. In order to do this, we make use of the 
low-energy theory (non-linear sigma model [3ll ]) for dis- 
ordered graphene with decoupled valleys derived in Ref. 
[iU (see also Ref. (32[)- The model is separated into two 
independent sectors corresponding to the two valleys. In 
each sector the action has the form 



S[Q] = \ Str 



.%x (VQ) 2 



9xy ± ^ )QV x QV a C 



(3) 

The field Q is the 4x4 supermatrix operating in Fermi- 
Bose and advanced-retarded spaces. The two parameters 
of the model, g xx and g xy , are longitudinal and Hall con- 
ductivities per one valley and per spin component mea- 
sured in units e 2 /h. The 'Str' operation implies the su- 
pertrace in all indices of the matrix along with real-space 
integration. This action differs from the usual sigma 
model in quantizing magnetic field [2CJ by the addition of 
±1/2 to g xy . This additional contribution arises due to 

Hil. It 



the quantum anomaly of Dirac fermions [25 
is the only reminiscent of the Dirac nature of excitations 
that survives at large scales and influences the critical 
properties. The signs in front of the anomalous terms 
1/2 are opposite for the two valleys. This ensures the 
global parity symmetry (x 1— > —x, K 1— > K') of the total 
action. 

The second term of the action has a topological 
nature: Sti(QV x QV y Q) = 8it:N[Q} with N[Q] taking 
only integer values. This term gives the imaginary part 
of the action ImS^Q] = 6N[Q] with the vacuum angle 
8 = 2ng xy ± n. 

The initial values of g xx and g xy are determined by 
the corresponding Drude expressions (see Supplementary 
Materials). The quantum corrections that establish lo- 
calization, and hence the QHE, are the result of renor- 
malization of the action Eq. (|3|). The renormalization 
flow of g xx and 8 was proposed in Ref. 2^, 34 1 . We plot 
this flow schematically in Fig. Q] by dotted lines. The 
effective theory (|3]) is invariant with respect to the vac- 
uum angle shift 8 1— > 8 + 2ir hence the flow pattern is 
periodic function of g xy . Transitions between quantum 
Hall plateaus occur when the value of 8 passes through 
an odd multiple of w. Owing to the anomalous contri- 
bution in Eq. ([3]), this is the case at zero filling factor 
when g xy — 0. Thus we have shown the validity of the 
odd quantization series Eq. J]} in the case when disorder 
does not mix the valleys. The absence of anomaly would 
have led to a plateau rather than the transition at n = 
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FIG. 1: Renormalization group flow of a xx and a xy in 
graphene with decoupled and mixed valleys. Dotted/dashed 
lines are separatrices of the flow for graphene with decou- 
pled/mixed valleys. Open circles are unstable fixed points cor- 
responding to quantum Hall transitions. Stable fixed points 
(plateaus) are shown as disks. Two solid curves demonstrate 
a possible flow towards even- and odd-plateau fixed point for 
a model with weakly mixed valleys. Each curve has a cusp 
when the running scale reaches Z m i x - 

similarly to ordinary QHE. 

Physically, the step of Hall conductivity between 
plateaus is due to a critical delocalized state which is 
exactly at the Fermi energy when 9 = ir. All other states 
are localized and do not contribute to either longitudinal 
or Hall conductivity. The value of longitudinal conduc- 
tivity exhibits a peak at the transition point with the 
maximum value 

ffxi = 4x 9 ^ 2e 2 /h, (4) 

where g^ is the longitudinal conductivity for the ordinary 
quantum Hall effect (known to be in the range g^ ~ 
0.5 -7- 0.6 from numerical simulations (35| ) and the factor 
4 again reflects the valley and spin degeneracy. Equation 
(|4|) agrees with the experimental value found in strong 
magnetic field at the Dirac point [4j, |5| |8( . 

Valley mixing: ordinary quantum Hall effect. — 
Let us now turn to the case when a weak valley mix- 
ing is present. For instance, charged impurities scatter 
electrons between valleys at some small rate as com- 
pared to the intra-valley scattering rate r _1 . The total 
action of the system will then be perturbed by the small 
coupling between matrices Qk and Qk' corresponding 
to the two valleys 

S[Q K , Qk'} = S[Q K ] + S{Q K >] + — Str Q K Q K ,, (5) 

7~m ix 

where S[Qk,K'] is given by Eq. ©, and p is the density of 
states at the Fermi level (see Supplementary Materials for 
the derivation). This perturbation is relevant and leads 
to the constraint Qk = Qk' in the infrared limit. The 
corresponding valley-mixing length is determined by the 
relation l m { X /l ~ {jmix/T) 1 / 2 . The ultraviolet scales I and 
t are given by the effective mean-free path and time; in 
strong magnetic field (for low-lying Landau levels with 
| AT | ~ 1), the length I is of the order of the magnetic 
length: I ~ Is = \fh~cj 'eB and the mean free time r. 
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FIG. 2: Quantum Hall effect in graphene with smooth disor- 
der at zero temperature. Hall conductivity as a function of 
the filling factor: odd (decoupled valleys, dashed line) vs nor- 
mal (weak valley mixing, solidline) quantization. Inset shows 
the energy dependence of the density of states. The state in 
the center of Landau level is delocalized (dashed lines) when 
the valleys are decoupled. The valley mixing splits this delo- 
calized state (solid lines). 

The index is z = 2 for non-interacting electrons (diffusion 
propagation) and in the case of short-range interaction 
[36j. A different value z emerges in the case of Coulomb 
interaction [33]; experiments [38| yield z ~ 1. 

At a scale larger than Z m ; x , we have Qk = Qk 1 and 
the topological terms with anomalous factors ±1/2 cancel 
in Eq. We end up with the unitary sigma model 

for the normal QHE [l9( with 9 — Aitg xy and ordinary 
quantization of Hall conductivity 

cr xy = k 2e 2 /h, keZ. (6) 

A delocalized state at the center of each Landau level 
is doubly degenerate when the valleys are decoupled. A 
weak valley mixing leads to a small splitting of the de- 
localized state within a single broadened Landau level. 
The new even plateau appears between the two odd ones 
when the chemical potential lies between the two split 
delocalized states (see Fig. [5]). The longitudinal conduc- 
tivity o~ X x has two separated peaks 2 x g'^ ~ e 2 /h in this 
case (here the factor 2 is due to the spin degeneracy). It 
is worth mentioning a similarity of the splitting of the 
anomalous QHE and the splitting of delocalized states 
by spin-orbit (spin-flip) scattering in a spin-degenerate 
ordinary QHE |39j |. 

The flow of o~ xx and o~ xy for both cases of decoupled 
and mixed valleys is shown in Fig. [1] For weakly mixed 
valleys (solid lines), a crossover occurs between these 
two flows at the length / m ; x . The even plateaus are 
much shorter than the odd ones ((T|) provided the val- 
ley mixing is weak. If the valleys arc completely decou- 
pled, the quantum Hall transition between two succes- 
sive odd plateaus has a finite width determined by the 
temperature-dependent dephasing length l v . The states 
close to the center of Landau level are localized at length 
that diverges as h oc ~ l(5n)~ v where Sn is the deviation 
of the filling factor n — 2itl 2 B n e from the transition point 
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FIG. 3: Quantum Hall transition at finite temperature. A 
double step in a xy and a double peak in a xx (solid lines) 
require low temperature, T < ft/r m i x . Otherwise a single 
broadened quantum Hall transition is seen (dashed lines). 

and v ~ 2.3 is the conventional quantum Hall critical in- 

l/u 



dex. The width of the transition is then 5r 



If the valley- mixing length Z m i X is larger than l v , the even 
plateaus will be totally smeared — the splitting between 
critical states is smaller than the delocalized energy re- 
gion around them. The even plateau becomes visible 
at sufficiently low T, when l v exceeds / m ; x (see Fig. [3]). 
Therefore, the width of this new plateau is 



(7) 

For Coulomb impurities, we estimate a typical value of 
splitting as Sn cvcn ~ 0.05 for the lowest Landau level and 
~ 0.1 for higher levels (see Supplementary Materials). In 
experiment, the temperature should be low enough in 
order to resolve the quantum Hall transition splitting, 
T < h/T miK . This implies T < 100 mK for the lowest 
Landau level and T < 1 K for higher levels. These val- 
ues are in reasonable agreement with weak localization 
measurements in low magnetic field ■ At higher tem- 
peratures, a broadened double step of Hall conductivity 
will be seen instead of two split transitions (Fig. [3]). 

Recent numerical studies [27[ demonstrated the split- 
ting of quantum Hall transition in graphene with a com- 
bination of potential and bond disorder. At the same 
time, the model with only potential disorder was found 
to show only the odd QHE in Ref. 27J . On the contrary, 
our consideration yields the existence of even plateaus in 
this case but with a narrower plateau at zero filling fac- 
tor (see Supplementary Materials). The zeroth plateau 
arises due to Landau level mixing which was discarded 
in Ref. 



27] 



Two other mechanisms, apart from intervalley scat- 
tering, can establish the even quantum Hall plateaus, 
Zeeman splitting and electron-electron interaction. Zee- 
man effect is weak in graphene; however, in Ref. Q 
the zero plateau that emerged in high magnetic field 
was attributed to this mechanism. An alternative - 
Stoner - mechanism was advocated in Ref. Q. Indeed, 
the repulsive interaction between electrons may result 



in the Stoner instability [4l|, |42| giving rise to spon- 
taneous breaking of spin and/or valley symmetry. Let 
us note that this instability would completely split the 
Landau level leading to the formation of even quantum 
Hall plateaus with the width comparable to that of odd 
plateaus as the magnetic field or electron mobility is in- 
creased 41] . This can be used to experimentally dis- 



tinguish the Stoner splitting from the disorder-induced 
splitting analyzed in the present work. 

Chiral disorder: "classical" QHE. — So far, we 
have considered the situation of a generic disorder within 
each valley. In Ref. [2l| it was shown that once the chiral 
symmetry Co is preserved by the disorder (e.g. ripples), 
the longitudinal conductivity at zero energy is exactly 
4e 2 /wh. External magnetic field also does not violate the 
chiral symmetry and hence does not change the value of 
conductivity [43[ . This leads us to the conclusion that the 
quantum Hall transition occurring at zero filling factor is 
modified by the presence of Co symmetry, since a xy — 
4e 2 /irh differs from the universal value ((4|) characteristic 
for a normal quantum Hall transition, a xx sw 2e 2 /h. 

A general form of chiral disorder in a single valley is 
a random (Abelian) vector potential A(r). The zeroth 
Landau level remains exactly degenerate in this situa- 



tion 



441 ] , as follows from the Atiyah-Singer theorem [45j . 
Moreover, one can find explicitly the wave functions at 
zero energy (see Supplementary Materials). The exact 
degeneracy of the Landau level implies the absence of 
localization. When the chemical potential lies at zero 
energy, the system behaves exactly as if it were clean. 
This means that the Hall effect is classical rather than 
quantum with a linear dependence of Hall conductivity 
on electron concentration n e 



&xy = n e ec/B = n 4e 2 jh. 



(8) 



This classical dependence holds for filling factor within 
zeroth Landau level, |n| < 1/2. The longitudinal conduc- 
tivity remains constant, a xx = 4e 2 /irh, in this case. The 
behavior of the Hall conductivity is shown in Fig. [U 

Let us now include a weak valley mixing maintaining 
the Co chiral symmetry. For instance, this is the case 
when the main disorder due to ripples is accompanied by 
rare dislocations. Let us recall that in the case of random 
scalar potential, the intervalley mixing leads to the split- 
ting of the quantum Hall transition into two with a small 
&xy — plateau in between. The longitudinal conductiv- 
ity is zero in this case. One could thus expect a similar 
behavior for chiral disorder. However, the conductivity 
at n = remains Ae 2 /irh according to the result of Ref. 
2ll as long as the chiral symmetry is preserved. This im- 
plies no QHE plateau. 

How does it happen that the valley mixing does not 
induce a quantum Hall plateau around the Dirac point? 
The answer is the same as for the Abelian random vec- 
tor potential discussed above: the zeroth Landau level 
remains exactly degenerate. The disorder we consider 
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FIG. 4: "Classical" QHE in graphene with chiral disorder 
(random vector potential). Chiral symmetry protects degen- 
eracy of the lowest Landau level (Inset: delta-function in the 
density of states). Hall conductivity is a linear function of 
carriers concentration while the lowest Landau level is be- 
ing filled. In Abelian case (ripples) only odd plateaus appear 
away from zero energy (dashed line). Non-abelian gauge dis- 
order (dislocations) split quantum Hall transitions as shown 
by solid line. 

corresponds to a random non- Abelian vector potential, 
A(r), which is a matrix in the valley space. The degen- 
eracy of the N = Landau level is a direct corollary 
of the Atiyah-Singer theorem 45| . An explicit construc- 
tion of zero-energy wave functions [4o| is given in Supple- 
mentary materials. The Hall conductivity again behaves 
classically within the zeroth Landau level, a xy — 4ne 2 /h 
for |n| < 1/2, but the other quantum Hall transitions, 
away from n — 0, split into pairs with narrow plateaus in 
between (see Fig. 2]), in the case of weakly mixed valleys. 

The observation of a narrow quantum Hall transition 
in graphene at n — seems to indicate that the dominant 
scattering mechanism is provided by long-range potential 
impurities rather than by ripples or dislocations. This is 
in agreement with the observed value of the zero-B mini- 
mal conductivity at the Dirac point which is appreciably 
larger than 4e 2 /irh expected for a random vector poten- 
tial. On the other hand, very recent experimental study 
of quantum Hall gaps in graphene [13] revealed that the 
lowest Landau level is significantly narrower than other 
Landau levels. This can be a signature of preserved chi- 
ral symmetry, suggesting that the main scattering mech- 
anism is due to ripples in the samples studied in Ref. 
471 ] . The quantum Hall measurement would provide a 
powerful test of this conjecture. 

Double-layer graphene. — Let us turn to the QHE 
in double-layer graphene. We limit our consideration to 
the case of disorder which does not mix the two valleys. 
The single-valley Hamiltonian of double-layer graphene 
reads 14811: 



H=±-[a x (pl 
2m 



P 2 y) + 2(TyP X Py]. 



(9) 
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FIG. 5: QHE in a double-layer graphene with smooth disorder 
(decoupled valleys) . Degeneracy of the lowest Landau level is 
twice larger than for other levels. Double step at zero filling 
factor (dashed line) is split when the disorder has finite corre- 
lation length d. Inset shows the density of states and positions 
of delocalized states (solid lines). Two such states within the 
lowest Landau level remain degenerate (dashed line) in the 
limit d — > oo. 

eB/mc. The two lowest levels, N = and N = 1, are de- 
generate. The corresponding wave functions are spinors 
in the sublattice space: (0, V'o) 7 ' and (O,-0i) T , respec- 
tively, where ipff is the wave function of N-th Landau 
level in a normal metal. 

In the presence of a generic disorder within each val- 
ley we have the same action Eq. © but with doubled 
couplings. The anomalous contribution to the topologi- 
cal term gives now 9 = 2it rather than it at zero energy. 
This implies complete localization and hence a plateau 
at n = 0. However, in experiments a plateau transi- 
tion with the double step in a xy at n — is observed 
instead. This can only happen if the disorder does not 
mix the two degenerate Landau levels with N = and 
N = 1. The only possible reason of the lack of mixing 
is the smoothness of disorder on the scale of magnetic 
length I g = y/Hc/eB. Indeed, the wave functions of the 
two Landau levels are orthogonal and concentrated in 
the area of order l B . If the disorder potential is almost 
constant in this small region, the corresponding matrix 
element is suppressed due to the orthogonality of wave 
functions. More specifically, assuming the disorder corre- 
lation length d^$> Ib, the mixing rate of the two Landau 
levels is found as: r^ 1 ~ r _1 (/s/(i) 2 . Comparison of tqi 
with the time needed for localization gives us the width 
of the zeroth plateau (see Fig. [5]) 



5no ~ (h/d) 



2/vz 



(10) 



The Landau levels are €n = hu> c ^ N{N — 1) with 
the conventional definition of cyclotron frequency lu c = 



To resolve this plateau, one should satisfy an upper 
bound on temperature, T < H/tqi- 

It is worth noting that the experimentally measured 
double step of the Hall conductivity at n — can not 
be automatically explained by the charged impurities in 
graphene. A random potential due to charged impurities 
has no characteristic length d. The only scale associated 
with such disorder is the screening length that is of the 
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order of Is (the stronger is the magnetic field, the larger is 
the density of states in Landau level, the more efficient is 
the screening) [i^] ■ The experimental observation of the 
double step in double-layer graphene thus suggests that 
an additional scale exists characterizing the smoothness 
of disorder. It might be caused by impurity correlations 
or, else, by their separation from the graphene layer. 

So far, we have considered the QHE in a single valley 
of a double-layer sample. If we include an intervalley 
scattering in our model, than the 4e 2 //i quantum Hall 
steps will further split, similarly to single-layers studied 
above. As a result, the conventional QHE with 2e 2 /h 
steps will be fully restored. 

Summary. — In this paper we have developed the 
theory of integer QHEs in graphene. The Landau level 
structure by itself is not sufficient to determine the form 
of the QHE. Anomalous QHEs in graphene are due to 
special character (symmetry) of disorder. In particular: 
(i) a smooth random (scalar) potential which does not 
couple the valleys gives rise to the odd QHE, Eq. ([1]) 
and dashed line in Fig. [2J (ii) the valley mixing splits 
the odd quantum Hall transitions and restores the or- 
dinary Hall quantization, Eq. ([6|) and solid line in Fig. 
[21 For weakly mixed valleys the crossover from the odd 
to ordinary QHEs occurs at parametrically low but still 
accessible temperatures, T < r~ ix ~ 100 mK; (iii) rip- 
ples or dislocations (random vector potential preserving 
the chiral symmetry) lead to a "classical" QHE, Eq. ([8|) 
and Fig. 2J around the half filling; (iv) in double-layers, 
a double-step QHE transition at n = arises for disor- 
der smooth on the scale of Ib- Experiments on QHE in 
graphene thus provide information about the nature of 
disorder. 
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Disordered graphene in strong magnetic field 

Here we present the calculation of the averaged Green 
function and of the density of states in disordered 
graphene in the presence of a strong external magnetic 
field. We will use the results of this calculation below for 
the derivation of the non-linear sigma model. We assume 
h = 1 from now on. 

Self-consistent Born approximation 

Let us start with the self-consistent Born approxima- 
tion (SCBA) approach. We assume single layer graphene 
with Gaussian (5-correlated disorder and consider first the 
simplest case of potential disorder characterized by a di- 
mensionless coupling constant a (it corresponds to a$ 
in Ref. 21]). This type of disorder does not produce any 
valley mixing, so that we can use the single- valley Hamil- 
tonian. Intervalley scattering processes will be included 
later. 

For the single- valley case, the Green function is a 2 x 2 
matrix in the sublattice space. In the presence of mag- 
netic field, the disorder-induced self-energy matrix has 
two distinct components, Si, 2, yielding the Green func- 
tion, 



G{e) 



ei 

V TC A 



V TT- 
(2 



tt± = p x ± ip y + - (A x ± i Ay 



where £1,2 = e — £1,2- 

The calculation of inverse matrix is straightforward in 
the basis of Landau levels. Owing to the fact that the 
disorder is (5-correlated, the SCBA equation involves only 
the Green function at coincident points. The latter is 
independent of a particular gauge and reads 



G(e;r,r) 



C2 E 

N' = l 







OO 



V 



ei E , 

N'=0/ 



The matrix SCBA equation S(e) = 2-Kv$aG(€\ r, r) de- 
termines two self energies Si ,2, 



aw 
~2 



2 I £2 E 

JV'=1 
00 

ei E 

N'=0. 



1 



eie 2 - u&N' 



(12) 



This equation was analyzed numerically in Ref. [131 ]. In 
the absence of magnetic field, the sum over Landau lev- 
els replaces by the integral and the well-known graphene 
SCBA equation is reproduced [Hoi ]. 



We are interested in the case of strong magnetic field, 
when Landau levels are well separated. Let us focus on 
a particular Ath level. Although SCBA gives the exact 
shape of the density of states only for N ^> 1, it yields a 
parametrically correct estimate for the height and width 
of the Landau level peak for all N . 

Consider first the case N 7^ 0. Then the largest term 
in the sums in Eq. (fl~2"|) is the one with N 1 — N. We 
estimate the sum of all other terms by replacing it with 
the corresponding integral. (It is worth noting that, con- 
trary to the case of a normal metal, the contribution of 
far Landau levels with N' 7^ N can not be neglected. In 
graphene, the density of states grows linearly with en- 
ergy; in strong magnetic field, this leads to Landau level 
separation decreasing as lu c / \fW' '. As a result, the contri- 
bution from high Landau levels to the self energy should 
be retained.) The difference between Si and S2, which 
originates from the term with N' = in Eq. (fT2|) . is im- 
material for N 0; we will use a unified notation S for 
them. We further simplify the equation by employing the 
inequality |e — ejv|, |S| <C and obtain 



S = 



4(e - e N - S) 



a(e-S)ln— . (13) 



The logarithmic divergence is cut by the graphene band 
width A. 

The effect of magnetic field is encoded in the u)\ term 



in Eq. (fl3|) . If this term were absent, the result would re- 
produce the well-known disorder-driven renormalization 
of the energy 21 1 



e = e - Re S = e/Z, Z=l-aln(A/e 



(14) 



It is instructive to express the solution of the full equation 
(I13p in terms of this renormalized energy e 



S= £+^ ± ^2 



CAT 



(15) 



The appeared parameter 7 determines the imaginary part 
of the self energy in the center of Landau level 



(16) 



2 A /l-aln(A/e JV ) 



The density of states within each Landau level has 
a standard form of semi-circle. The identity p(e) — 
— 7T _1 Imtr G R (e; r, r) together with the self-consistency 



9 



equation yields 



Im 



Si 



2ir 2 v 2 



2iv 2 v 2 a 
_ V4f* - ( e - e N Zy 
Att 2 1 2 b j 2 



(17) 



The last expression contains two parameters: the renor- 
malization factor Z determines the rescaling of Landau 
levels according to Eq. (fT4|) and the electrons scattering 
rate 7 = jZ gives the Landau level width. 

The result (fT7|) together with Eq. (fTBj) provides the 
following criterion for the separation of Landau lev- 
els within SCBA: the Nth level becomes isolated when 
u c Z/y/N > 7. The stronger is the magnetic field the 
larger is the number of isolated Landau levels in the vicin- 
ity of the Dirac point. 

The solution of the SCBA equation is qualitatively dif- 
ferent for the lowest Landau level, N = 0. The distinc- 
tion between Si and E2 is now crucial, since the term 
with N' = is absent in the sum (fT2|) for Si. Replacing 
the sums over non-zero levels with the integrals, we get 



Si = -a(e - E 2 ) ln(A/w c ), 

9 

£ 2 = 



2(e-E 2 ) 



a(e — Ei) ln(A/w c 



(18) 



We express the solution of these equations in terms of 
renormalized energy e according to Eq. (114|) 



E 2 = - ± i 



Ei = — aln(A/w c ) 



€ -±i 
2 



(19) 



We denote the imaginary parts of E12 at the center of 
Landau level by 71 2 



a ln(A/w c 
1 



^/2[l-a 2 ln 2 (AK)] 



(20) 



The electron scattering rates for the two sublattices are 
given by 71 _ 2 = ■^71,2 with Z from Eq. (| 14[) . The rate 
72 has the same order of magnitude as for a non-zero 
Landau level (fT6|) while 71 is somewhat smaller. This 
is a manifestation of the fact that the lowest Landau 
level wave function has its support in the sublattice B. 
With the opposite orientation of magnetic field, the wave 
function will be in A sublattice and 71 > 72. In the 
second valley the situation is reversed. 

Calculating the density of states at the lowest Landau 
level separately for two sublattices, we find 



For both sublattices, the width of the zeroth Landau level 
is 7 = 72Z, which determines the width of the total den- 
sity of states p = pi + p^- 

Ballistic renormalization group 

As we discuss in the previous section, high Landau 
levels produce logarithmic corrections to the low-energy 
properties of the system. The SCBA takes these correc- 
tions into account only partially. The systematic way 
for summing up such logarithms is the renormalization 
group (RG) formalism 0, [H, 01 . Below we develop 
this approach for the case of strong magnetic field. As 
we demonstrate, the results are qualitatively similar to 
the SCBA, but differ quantitatively. 

In the simplest case of diagonal Gaussian disorder a, 
the starting point for the renormalization group is the 
fcrmionic action 



d 2 r [-itp(e + iOA - H)ip + nv 2 ,a('ip'ip) 1 



(22) 

Here H is the single- valley Dirac Hamiltonian. The field 
-0 is an 8-supervector with the structure in the inner 
AB space (sublattices) of the Hamiltonian H, retarded- 
advanced (RA) space, and Bose-Fermi (BF) superspace 



22j. We use standard notation A = diagjl, — 1}ra- 
The doubling of variables in the RA space is needed for 
the calculation of averages involving both retarded and 
advanced Green functions, e.g. conductivity. In the bal- 
listic regime that we consider here, the distinction be- 
tween retarded and advanced propagator is immaterial. 

The renormalization procedure eliminates fast degrees 
of freedom thus reducing the cutoff energy A — * A/L. 
The parameters of the action (|2"2"|) are then rescaled ac- 
cording to [21 1 



da 
dhxL 



= 2a 2 , 



de 
d\nL 



(23) 



To study the properties of N ^ Landau level, we stop 
the renormalization at L = A/ejv, when the running cut- 
off reaches the observation energy. The new parameters 
are 



a 

Z2> 



e 



l-2aln— . (24) 



Now we employ the SCBA equation (fl3|) with the renor- 
malized parameters. The logarithmic term is absent as 
long as the running cutoff equals e^r after renormaliza- 
tion. The SCBA equation involves a single Landau level 
and yields the renormalized self energy E(e). In the cen- 
ter of Landau level, the imaginary part of E is 



7 



a/2. 



(25) 



pi\ _ fcdn(AK)j V 4V-' 



4:TT 2 v 2 aZ 



(21) 



The energy dependence of Im E gives the renormalized 
density of states /5(e) . In order to calculate the observable 
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density of states, we use the identity p/p = de/de and 
obtain 



y/aus 2 - (e - e N Z) 2 
2ir 2 vZa 



(26) 



This result has the same form as the result of SCBA (fT7|). 
but the parameters 7 = 7Z and Z are modified. 

At the lowest Landau level, we use equations (fl"8]) 
with renormalized parameters and omitted logarithmic 
terms. We calculate two self-energies, Si. 2; their imagi- 
nary parts at e = are 



7i =0, 



72 



(27) 



The renormalized density of states is concentrated in the 
sublattice B, 



Pi = 0, 



(>2 



\/2&uj. 



(28) 



In order to find the observable densities, we have to mod- 
ify our RG scheme. Different values of /5 12 call for intro- 
ducing the two different energies, ei.2, in two sublattices. 
The equations for these energies have the form 



ae 2 - 



dlogL 

The solution reads 

1 r . 
£1,2 = - [ei,2{Z + Z 



dlogL 



The connection between p and p has the form 



11=1,2 



Pp.- 



(29) 



(30) 



(31) 



Using Eqs. ([25]). (j3T)]) . and ([3T]), we obtain the resulting 
density of states at the lowest Landau level in two sub- 
lattices 



aln(A/w c ) ] y/2au>*-e 2 
-aln(A/w c ) J 



A-k 2 v 2 ol 



(32) 



This result is also similar to its SCBA counterpart Eq. 
(f2"Tj) with slightly modified parameters. The width of 
the Landau level is 7 = %Z . The lowest Landau level 
becomes isolated when to c exceeds Ae _1 / 2a . 

We can further improve the result by employing the 
exact density of states at the lowest Landau level found 
by Wegner [49(. After eliminating all non-zero Landau 
levels with the help of RG, we find ourselves in the situa- 
tion when the approach of Ref . [Zi| is directly applicable 
and yields 



92 



F{1/ 12) 
2^ 2 /| 7 2 ' 



F(x) 



l+[(2/^)f X ey 2 dy}- 



(33) 



By substituting this result into Eq. ([31]) , we calculate the 
observable density of states 



a ln(A/w c ) 
1 - aln(A/w c ) 



(34) 



This improves the result Eq. (f3"2"| by replacing the semi- 
circle function F(x) = \Jl — x 2 /4 with the exact lowest 
Landau level shape Eq. ([33]). 

Charged impurities 

So far, we have considered finite-range disorder. The 
model of long-range charged impurities can also be 
treated in the framework of SCBA once the screening 
is taken into account. For low- lying Landau levels, the 
screening length is of the order of magnetic length, which 
is the only scale in magnetic field [42j. For high Lan- 
dau levels, the screening occurs at a scale of the electron 
wavelength. If the dimensionless parameter characteriz- 
ing the interaction strength is small, r s — e 2 /v Q x *C 1 (x 
is the dielectric constant), the screening can be control- 
lably treated within the random phase approximation. In 
a more realistic situation, r s ~ 1, the results for charged 
impurities are valid up to a numerical factor of order 
unity in the definition of the effective disorder strength 



l/N, N ^0, 
1, 



N = 0. 



(35) 



The density of states follows from the SCBA equations 
{13} or (pi]) with a from Eq. (J35j) and without the log- 
arithmic terms [i.e., ln(A/e) is replaced by a number of 
order unity]. The result has the form of Eq. |26|) with 
Z = 1. The absence of the logarithmic terms is due to 
the suppression of scattering off Coulomb impurities at 
large momentum transfer (that is, transitions involving 
far Landau levels are ineffective) . The lack of hard scat- 
tering also leaves no room for ballistic rcnormalization. 

Derivation of the sigma model 

Sigma model in a single valley 

Non-linear sigma model is an effective low-energy the- 
ory describing soft modes of the system, diffusons and 
Cooperons [31]. In the absence of valley mixing, the 
sigma model for graphene in zero magnetic field was de- 
rived in Ref. [22I ]. Here we generalize this derivation, al- 
lowing for the magnetic field within a single valley. Then 
we will also include the intervalley scattering. 

We start the derivation from the fermionic action (|22D . 
The RA structure of the fields will play a crucial role 
in the sigma model. Our calculation is based on the 
SCBA approach outlined above. The more rigorous RG 
calculation can also be used (as in the zero-B case [Til ]) 
as a basis for the sigma model, leading to the same form 
of the theory. 
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The ('ipip) 2 term in Eq. (f22|) is decoupled with the help 
of an auxiliary 8x8 supermatrix field R. Subsequent 
Gaussian integration over ip yields an effective action in 
terms of R, 



S[R] 



StrR 2 



Strln [e-H-R]. 



(36) 



The soft modes of the system, that sigma model deals 
with, describe the fluctuation near the saddle point of 
S[R\. This saddle point is determined by the self- 
consistency equations (fT^j) with the self energy E replaced 
by the matrix R. We separate the real and imaginary 
parts of the self energy 



R = Re E + iTA, 



(37) 



where T is the matrix of renormalized scattering rates, 
f = diag{ 71,72} ab, given by Eqs. (fl6|) or ([20]) . A whole 
saddle manifold can be generated from the solution ((37|) 
by a uniform rotation T that commute with the Hamil- 
tonian H. As a result, the matrix A in the imaginary 
part of Eq. ([37]) replaces with Q = T^AT. The 4x4 
matrix Q operates in RA and BF spaces and obeys the 
constraints strQ = and Q 2 = 1. We rewrite the action 
(|36|) in terms of Q omitting the first term which produces 
an unphysical constant, 



S[Q] = Strln e-H + iTQ 



(38) 



The real part of the self energy is included in e which be- 
comes an AB matrix. Effective low-energy action (sigma 
model) is a result of the gradient expansion of Eq. ([38]) . 
This expansion is a non-trivial procedure in view of the 
topology of the saddle manifold [2(j. Furthermore, the 
Dirac nature of electrons in graphene gives rise to extra 
anomalous contributions to the sigma-model action [22| . 

The approach of Ref. [22j is directly applicable to the 
derivation of the sigma model in magnetic field. The key 
feature of this approach is a special form of boundary 
conditions involving the mass term, 77103, m the Hamil- 
tonian. Assuming the mass is zero in the bulk of the 
sample and gradually increases up to some large value 
M near the boundary, we get the sigma-model action 

M 



S[Q] = i Str 



-^(VQ) 2 + — QV x QV y ( 
Z Ztt 



(39) 
(40) 



with the topological angle 

9 = 2n[g I xy (0)+g I x I y (0)-g I x I y (M)}. 

The parameters of the model are determined by the stan- 
dard Kubo expressions 

9** = ~~ Tr [j x (G R - G A ) Jx (G R - G A )] , (41) 

<4 = -\ Tr [MG R G A ) Jy (G R + G A )] , (42) 

(43) 



J:ry 



= -T T [( XJy -yj x )(G R -G A )] 



Trace in the last equation is divergent and requires a reg- 
ularization. This happens because g xy accounts for the 
contribution of edge modes to Hall conductivity. That is 
why we have to specify boundary conditions in order to 
find g". 

The dependence of g xy on boundary conditions shows 
that the very notion of the single- valley Hall conductiv- 
ity can not be properly defined. The observable Hall 
conductivity 



9xy 



9xy 



2V xy 



9xy\ B-t-1 



(44) 



always includes contributions from both, mutually time- 
reversed, valleys implying a cancellation of divergences 
in Eq. (|4"3"|). Considering the Hall conductivity per val- 
ley, one usually means a half of the total observable Hall 
conductivity. This corresponds to a certain regulariza- 
tion requiring g xy = at the Dirac point. 

At the same time, the value of 9 in the sigma-model 
action is well-defined (modulus 2ir) even within a single 
valley as long as contains a difference of two g xy quan- 
tities (|40l) . At the boundary, the introduced mass M is 
large, so we can neglect energy and magnetic field there 
and obtain fl^(M) - c^(0) = (1/2) sign M. This pro- 
vides the anomalous topological term in the sigma model 
Eq. with 8 = g xy +l/2 signM. The sign of the anoma- 
lous term (sign M here) is immaterial as it only changes 
the action by an integer multiple of 2ni. 

It is worth emphasizing that the localization or critical- 
ity is the property of the bulk theory and does not depend 
on the boundary condition. Nevertheless, similarly to the 
ordinary QHE, introducing the boundary turns out to be 
a convenient way of deriving the field theory, since the 
action contains a topological term. The resulting the- 
ory however does not depend on whether a system with 
boundary or without it (say, on a sphere) is considered 
and on the way the boundary is implemented. Indeed, 
the final form of the topological term in Eq. ([3]) is repre- 
sented as a 2D integral over the bulk. Thus the bound- 
ary only facilitates revealing and exploring the intrinsic 
topological properties of Dirac fermions in the bulk of 
graphene. 

An alternative derivation of the sigma model for Dirac 
fermions employs non-Abelian bosonization [24j that 
does not require an introduction of boundary conditions. 
In bosonic language, disorder leads to constraint on the 
boson field reducing the chiral gauge symmetry group 
down to sigma-model manifold. The Wess-Zumino term 
in the bosonized action transforms into the anomalous 
topological term of the sigma model. This method was 
used in Ref. [l|| for graphene with mixed valleys. 

Intervalley scattering 

Let us now add an intervalley scattering term to the 
fermionic action (j2"2")l . The intervalley scattering due to 
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time-reversal invariant disorder is described by two cou- 
pling constants, /3± and [3 Z (see Ref. [2lj for details), 



X Str|/3_L [(tp^AK^^AK' + {ip^BKilp^BK' 
+ (3 Z [(tjj^AK^^BK' + {'H^)bk('H^)aK']\- 



(45) 



We will treat this term perturbatively within the SCBA 
scheme. This is equivalent to replacing a pair of ip fields 
with the corresponding Green function which, on the 
saddle-point level, is equal to the matrix Q: 



(#) 



K,K' 



^K,K'Qk,K 

2ixvl\a 



(46) 



The imaginary part of self energy is different in two val- 
leys, f K = diag{7i,7 2 }AB, T K ' = diag{7 2 , 7i}ab- After 
the substitution (146|) . the valley-mixing action acquires 
the form of Eq. ((SJ). We calculate the mean free time 
from the width of Landau level, r = 1/(47), an d obtain 
in the level's center 



2uj c ot 



27r/3 ± 7 1 7 2 +7r/3 z (7 1 2 + 7 2 2 ) 



(47) 



The form of the 5 m i x term is universal and does not 
rely on the particular disorder model. At the same time, 
the mixing rate r m i X is determined by microscopic non- 
universal mechanisms and depends on the disorder type. 
A potential disorder provides only the intervalley cou- 
pling P±. Using the SCBA results ([16j) and 1(20]). we find 



■ 4a [1 - aln(A/ejv)] 

2[l-a 2 ln 2 (A/w c )] 
tt/3 ± ln(A/w c ) : 



JV/0, 



N = 0. 



(48) 



In order to apply the ballistic RG approach, we first 
renormalize the action (|22[) including the mixing term 
()45|) . Assuming the inequality a ^> 0± z , we employ the 
simplified version of RG equations [l8j, \2l\ 



d\nL 



Aa(3 z , 



d(3 z 
d\nL 



= -2aj3 z + 2a/3 x (49) 



in addition to Eq. (|23[) . The renormalizcd couplings are 
then substituted into Eq. (|47[) with the parameters 7^2 
given by Eqs. ([25]) and (|27p . In terms of bare couplings 
a and (3±_ (potential disorder), the mixing time is 



4a 



7T (3 ± In A [1 _ 2a In A + 4*1 ]^ A] 



A^O, 



N = 0. 



(50) 

For the lowest Landau level, the RG rate 71 appears to 
be zero 27J since the wave function resides solely in the 



sublattice B. The valley mixing occurs only due to /3 Z 
disorder. For potential impurities, this coupling has been 
absent in the ultraviolet limit but is generated by the RG 
flow dHJ. 

The SCBA and RG results @SJ| and coincide up to 
a numerical factor of order unity once the Landau levels 
are well separated, i.e. in the range of our interest. The 
criterion of level separation is provided by RG calcula- 
tion: uj c > Ae^ 1 / 2 ". 

Estimation of T m i x 

Now we apply the model of screened Coulomb impu- 
rities to find the values of r and r mix . This model is 
most relevant for graphene experiments because it con- 
forms with both linear dependence of conductivity on the 
concentration of electrons J42j and with minimal conduc- 
tivity at the Dirac point [9j. 

The Born parameter of screened Coulomb impurities 
is given by Eq. (|35|) . Intervalley scattering involves large 
momentum transfer. Thus we neglect screening and es- 
timate (i±_ ~ WimpO 2 where a is the lattice constant. The 



valley mixing rate follows from any of Eqs. (|48|) or ([50] 
with logarithmic factors replaced by some numbers of or- 
der unity. 



l 2 B /a 2 N, N 0, 
(ntapo 2 )- 1 , A = 0. 



(51) 



Taking a typical magnetic field value of 20 T, we get a 
10% splitting of quantum Hall transitions for non-zero 
Landau levels. Another splitting of order 5% appears at 



zero Landau level if we estimate 



''imp 



4 x 10 n cm- 2 



from mobility measurements away from the Dirac point 

a. 

To observe the quantum Hall transition splitting, the 
temperature should be small enough. Namely, one should 
have T < l/r m i X . For non-zero Landau level at 20 T this 
gives an upper bound of ~ IK. The splitting of zero level 
becomes visible at smaller temperatures ~ 100 mK. 

Wave functions in chiral disorder 

Abelian case 

Abelian chiral disorder has a form of random vector 
potential. After a proper gauge transformation, any two- 
dimensional vector potential can be expressed as a curl 
of a scalar field <h(r) 



A x = -V, 



Ay V :t 



(52) 



This field 4> is uniquely determined by the magnetic field 
B(r) penetrating the system, V 2 = —B. Assume that 
the uniform part of magnetic field Bo, the one that es- 
tablishes Landau levels, is pointing up, Bq > 0. Then 
the function <f> grows at infinity as <j> ~ B^r 2 and all the 
zero-energy wave functions lie entirely in the sublattice 
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B. A possible set of such functions (up to a normalization 
factor) is 

*£(a:, y) = (x - iy) m exp(-e^/c). (53) 



where A± = A x ± iA y , d± — W x ± iV y , and g is an 
appropriate 2x2 matrix in the valleys space. The wave 
functions of the zeroth Landau level again lie in the sub- 
lattice B and have the form 



Non-Abclian case 

Non-Abelian vector potential has a matrix structure in 
the valleys space. An explicit construction of zero-energy 
wave functions is almost the same as above (46| : express 
the vector potential in the form 

A+ = - g- l d + g, A_ = gd^g' 1 , (54) 



*Z(x,y) = (x-iy) m gi , 2 (55) 



with <7i,2 being any of the two columns of the matrix g. 



